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In this work, we use three different numerical techniques to study the charge transport properties of a system in 
the two-level SU(2) (2LSU2) regime, obtained from an SU(4) model Hamiltonian by introducing orbital mixing 
of the degenerate orbitals via coupling to the leads. SU(4) Kondo physics has been experimentally observed, and 
studied in detail, in Carbon Nanotube Quantum Dots. Adopting a two molecular orbital basis, the Hamiltonian 
is recast into a form where one of the molecular orbitals decouples from the charge reservoir, although still 
interacting capacitively with the other molecular orbital. This basis transformation explains in a clear way how 
the charge transport in this system turns from double- to single-channel when it transitions from the SU(4) to the 
2LSU2 regime. The charge occupancy of these molecular orbitals displays gate-potential-dependent occupancy 
oscillations that arise from a competition between the Kondo and Intermediate Valence states. The determination 
of whether the Kondo or the Intermediate Valence state is more favorable, for a specific value of gate potential, 
is assessed by the definition of an energy scale To, which is calculated through DMRG. We speculate that the 
calculation of To may provide experimentalists with a useful tool to analyze correlated charge transport in many 
other systems. For that, a current work is underway to improve the numerical accuracy of its DRMG calculation 
and explore different definitions. 

PACS numbers: 73.63.Kv, 72.15.Qm, 73.23.Hk 



I. INTRODUCTION 

Interest in charge transport properties of carbon nanotubes 
(CNTjD has grown consistently in the last few years, in part 
because of their possible applications as single electron tran- 
sistors (SETs), 2 as well as their possible use in molecular 
computing. 3 -' A CNT SET can be manufactured by coupling 
a CNT quantum dot (QD) to metal leads 2 . In this case, be- 
low a certain characteristic temperature TV, the Kondo effect 4 
was first observed in CNT QDs in 2000. 5 Further investiga- 
tion showed that the main difference between the Kondo ef- 
fect in CNT QDs, and the same effect in other nano-structures, 
comes from the degeneracy of the orbital (chiral) states of the 
CNT QD. 6 Effectively, while in the one hand the Kondo ef- 
fect in, e.g., a lateral semiconducting QD is associated to the 
screening of a localized magnetic moment (in general, with 
SU(2) symmetry) by the spin density of the conduction elec- 
trons, in the other hand, in a CNT QD the two degenerate 
orbital states with different chirality (usually referred to as a 
pseudospin), together with the intrinsic spin degree of free- 
dom, give rise to the so-called SU(4) Kondo effect. This ex- 
otic state has been recently studied both experimentally^ and 
theoretic allyPKS Some of the many particular characteristics 
of the SU(4) Kondo effect are (i) a larger Kondo tempera- 
ture than in the SU(2) case, (ii) Kondo effects at two different 
fillings, namely, quarter-filling and half-filling, with diverse 
properties,^! (hi) a rich behavior under magnetic field,^ (iv) 
perfect entanglement of spin and orbital degrees of freedom, 6 
among others. 

Another interesting regime in the SU(4) Kondo Hamilto- 



nian can be obtained by the gradual 'mixing' of orbital states 
at the tunneling barriers, through tunneling that effects a flip 
in the pseudospin of the electron [henceforth, pseudospin con- 
serving tunneling occurs through a matrix element t', and non- 
conserving through t", see Fig. 1(a)]. When this mixing is 
total (i.e., t" /t 1 = 1), one reaches the so-called two-level 
SU(2) (2LSU2) state.^One constraint associated to the use 
of a CNT QD to observe SU(4) Kondo physics is the lack of 
control over the parameters of the system (intra- and inter- 
orbital Coulomb repulsions U and U' , respectively, and cou- 
plings to leads), making it impossible to freely navigate from 
the SU(4) to the 2LSU2 state. In contrast, recent advances in 
the lithography of semiconducting lateral QDs have allowed 
greater control over capacitively coupled double quantum dot 
(DQD) systems connected to independent leads. 12 In reality, it 
was in a DQD system that the first claim of observation of the 
SU(4) Kondo effect was made. 13 As long as progress in this 
field continues to be made, it is not unlikely that the 2LSU2 
regime will be accessible to experimental probing in the near 
future. Indeed, in the same way that the SU(4) regime has 
been shown to be quite robust regarding a difference between 
U and U', and some mixing of the orbital states (although 
strict SU(4) symmetry only occurs at one point in parame- 
ter space), 8 -* most of the relevant phenomena described here 
for the 2LSU2 state are shown to survive a deviation of if' /if 
from 1. It is for that reason that we elect to use the terminol- 
ogy SU(4) and 2LSU2 regimes, since as far as experiments 
go, their characteristic properties do not seem to be restricted 
to the corresponding point in parameter space (we choose to 
refer to these specific points as SU(4) and 2LSU2 states). 
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In this paper, using three different numerical techniques, 
namely the Numerical Renormalization Group (NRG)p^ the 
Density Matrix Renormalization Group (DMRG), 15 and the 
Logarithmically Discretized Embedded Cluster Approxima- 
tion (LDEC A)J-^I we analyze in detail the opposite e nd of the 
SU(4) regime, the above mentioned 2LSU2 regimel^ESl Us- 
ing the notation defined above, the difference between the 
SU(4) and 2LSU2 regimes resides in the value of the ratio 
v = t"/t': v = corresponds to the SU(4) regime (no 
orbital mixing) and v = 1 to the 2LSU2 regime (total or- 
bital mixing). The motivations for the study presented here 
are four-fold: First, in previous work by two of the current 
authors, a charge discontinuity and a transition from a 2Gq 
to Go maximum conductance was obs erved in the 2LSU2 
regime (Go is the conductance quantum). 10 17 Second, in these 
same works, a fine-structure of conductance peaks (2Go high 
and positioned at and around the particle-hole (p-h) symmet- 
ric point) was observed for parameter values in the 2LSU2 
regime, superimposed over a flat background of height Go- 
Third, a slave boson analysis of the local density of states 
(LDOS) showed a discontinuous change of the Kondo peak at 
the 2LSU2 regime. 9 All these three points mentioned above 
were never fully explained, some of them being thought to 
be related to convergence problems in the numerical methods 
usedfSl Here, we will show that this is not the case, and that 
they are associated to a band-decoupling at the 2LSU2 state 
and to an interesting competition between two distinct many- 
body ground states. Finally, the above mentioned advances 
in manufacturing of capacitively coupled DQDs may open the 
doors for the experimental analysis of the regime theoretically 
discussed here. Note that to facilitate the connection and com- 
parison with previous work by some of the authors, we will 
keep referring to the system under study as a CNT QD, al- 
though the regime of parameters to be analyzed ( v < 1) may 
not be accessible to CNT QDs. 

Another interesting aspect of the results presented here is 
that the points mentioned above (charge discontinuity, change 
from double- to single-channel charge transport, fine structure 
of the conductance, and discontinuity in the LDOS Kondo 
peak) can be more easily understood after a change of ba- 
sis from 'atomic' to 'molecular' orbitals (AOs to MOs) is ap- 
plied. Indeed, QDs have been dubbed 'artificial atoms', and 
systems of two or more coupled QDs have been dubbed 'artifi- 
cial molecules'. When the phenomena being analyzed depend 
more on the properties of quantum states involving more than 
one QD (a 'molecular orbital'), a change of basis from local 
states (individual QDs) to extended states, involving two or 
more QDs, should be done to help understand the phenom- 
ena being observed.^ In this work, the authors will present a 
situation where this change of basis is crucial for the under- 
standing of the numerical results. 

This work is divided as follows: Section II briefly presents 
a description of the model Hamiltonian, and the regime of pa- 
rameters to be analyzed. In addition, a brief description of the 
band decoupling occurring at the 2LSU2 state is described by 
introducing molecular orbitals (a detailed description is pro- 
vided in Appendix A). Section III presents charge vs. gate 
voltage results obtained with DMRG and NRG, at zero tem- 



perature. Charge oscillations in the MOs, as a function of gate 
voltage, are in perfect agreement between the two techniques. 
The NRG results are extended to finite temperatures, deter- 
mining the energy scale of the phenomena relevant for the 
^-dependent charge oscillations. In addition, to facilitate the 
discussion of the results, the idea of an effective gate poten- 
tial V* is introduced. Section IV presents LDECA results for 
charge, conductance, and spin- and charge-fluctuations, away 
from the 2LSU2 point, in good qualitative agreement with the 
DMRG results. Section V presents a discussion of the origin 
of the Vg-dependent charge oscillations observed in the MOs, 
and, more importantly, an explanation of its dependence on 
the coupling of the CNT QD to the leads. These charge os- 
cillations in the MOs reveal a subtle competition between two 
many-body states, namely, the Kondo state and the Intermedi- 
ate Valence state. A simple way of quantifying this competi- 
tion is obtained by using DMRG simulations to calculate the 
energy Tq of the system for a wide window of gate potential 
V g (see Eq. (12) for its definition and discussion thereafter). 
In addition, the DMRG calculation of To for the single impu- 
rity Anderson Model is compared to the Kondo temperature 
expression derived by Haldane. 20 A finite size scaling of the 
DMRG results shows that Tq tends to Haldane's expression in 
a narrow window around the p-h symmetric point, leading the 
authors to speculate that the definition and numerical calcula- 
tion of T may be helpful to study correlated charge transport 
for a broad region of gate potential, including the Kondo and 
Intermediate Valence regimes. Section VI presents our con- 
clusions. Appendix A presents details of the basis transforma- 
tion that results in the band decoupling essential to understand 
the properties of the 2LSU2 regime, and Appendix B shows 
results for the non-interacting model (U = U' = 0.0). 



II. HAMILTONIAN AND MOLECULAR ORBITALS 

A. Model Hamiltonian 

The SU(4) Hamiltonian used to model a CNT QD (or a ca- 
pacitively coupled DQDj-SI coupled to metallic leads is given 
by 



H to t — -ffmb + -ffband + -Hhyb, 

u 



H m b — 

<t; A— a,/3 



-n\„n\ s + V g nxo 



A— a,/3 i=l,a 
4_ iff C A -i-l ff ) +H.C, 



(1) 



(2) 



hyb 



A, A'— a,/3,*,tr 



d\ a c\>±i a +H.C. 



(3) 
(4) 



Our unit of energy is t, and all results were obtained for 
U = U' = t. Note that this Hamiltonian is p-h symmetric 
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FIG. 1: (a) Schematics of the full system studied in this work, rep- 
resenting either a CNT (in which case a and /3 stand for the orbital 
(chiral) levels of the CNT), or a capacitively coupled DQD (in which 
case a and j3 stand for two QDs that are coupled to each other only 
through a long range Coulomb interaction). In either case, the orbital 
levels of the QDs are each coupled to its own independent charge 
reservoir, through hopping t! [(green) line]. Coupling to a reser- 
voir with opposite quantum number occurs through hopping t" [(red) 
line]. The SU(4) state is achieved for t" = 0, and the 2LSU2 state is 
achieved for t" = t' . (b) Schematics of the usual left-right symme- 
try transformation that, in this case, decouples two bands from the 
interacting levels, (c) Final system, obtained after the second sym- 
metry transformation (see text). Now, the MOs |+) and |— ) are cou- 
pled to bands through hoppings t+ and i_ that depend on the ratio 
v — t" /t' . For v = 1 (2LSU2 state), level j— ) decouples from the 
electron reservoir, staying coupled to the rest of the system through 
its many-body capacitive coupling to level | +) . 



for Vg/U = —1.5. if m b contains all the many-body terms 
in the system; more specifically, A = a, f3 are two degen- 
erate orbitals associated to the wrapping mode (clockwise or 
anticlockwise) of the electron propagation along the axial di- 
rection of the CNT, while d aa (dp a ) annihilates an electron 
with spin a in the a (J3) orbital in the CNT QD and c\i a an- 
nihilates an electron with spin a in the i-th site of the A = a 
or f3 channel for i > or i < [right or left leads, respec- 
tively, see Fig. 1(a)]. Therefore, each lead has two bands, 
also named a and j3, after the notation used for the orbital 
levels in the QD. Justifications for this assumption are given 
elsewhere It is important to stress that the CNT QD is cou- 
pled to two independent metallic leads on each side (at left and 
right), which have the same chiral states as the CNT QD (and 
therefore the same quantum numbers). Equations (3) and (4) 
describe the leads and their coupling to the QD, respectively. 
More specifically, t attx — tp,p — t' and t a> p — tp ltx = t". 
To explore the 2LSU2 regime, we wish to perform numeri- 
cal calculations for v < 1. Following Ref. 9 we will de- 
fine a ratio <j> = tan^ 1 (z/), with which the hoppings between 
the QD and leads can be rewritten as t" = to sin(</>) and 
t! = to cos (</>), where to is kept constant, and t' and t" are 
indicated in Fig. 1(a). Using this definition ensures that a con- 
stant coupling T = 2irp (Ep) to (where p 1S the reservoir's 
DOS at the Fermi energy) between the QD and the charge 




FIG. 2: (a) DMRG results for charge variation against gate poten- 
tial: n+(V g ), (green) x symbols, n_(Vg), (blue) + symbols, and 
n a + up, (red) solid line (note that n a — np). The charging of the 
AOs is smooth, while the charging of the MOs has clear disconti- 
nuities. The parameters used are <j> = 7I "/4, to = 0.106, (7 = 1, 
and r = 0.045 (note that the particle-hole symmetric point of this 
model is at V g — —3(7 /2, in addition, the charge is plotted per spin 
species), (b) Comparison between NRG and DMRG n± vs. V g re- 
sults for the 2LSU2 state (4> = it/ 4, to = 0.025 (T = 0.0025), and 
U = 1). NRG: n+ [circles (black)], n_ [squares (red)]. DMRG: 
n+ [x (green)], n_ [+ (blue)]. Note the very good agreement be- 
tween the two techniques. DMRG results obtained with 40 sites and 
maintaining 200 states. NRG results obtained at the last iteration 
(N=80), and a discretization parameter A = 2.5 was used, keep- 
ing N a — 1500 states in each iteration, which allows a larger N„ 
to account for degeneracies (all the other NRG calculations in this 
work were done similarly). The small differences occur because the 
DMRG results are not totally converged yet. Note also [as in panel 
(a)] the interesting Vg -dependent charge oscillations involving the 
two MOs simultaneously, given that they are capacitively coupled. 
See text for definition of V gc and V g d- 



reservoirs is kept constant, while the ratio v varies. Note that 
v = it" = and t' = t Q ), the pure SU(4) state, implies 
= 0, while v = 1.0 (f = t" = t /V2), the pure 2LSU2 
state, implies cf> — it/ 4. In what follows we will assume, for 
simplicity, that U = U'f^ 

Figure 1(a) shows a schematic representation of H to t, to- 
gether with diagrams indicating the gradual decoupling of 
bands from the QD [Figs. 1(b) and 1(c)] as the symmetries 
of the system are exploited, finalizing with the construction of 
the MOs |+) and |— ). A brief description of these transforma- 
tions is given below and full details are provided in Appendix 
A. 



B. Molecular orbitals: band decoupling 

Given the left-right symmetry of the system, two 
bands decouple from the QDs, as indicated in Fig. 1(b), 
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through a symmetric-antisymmetric combination of creation- 
annihilation operators from left and right leads (without mix- 
ing a and (3 states). The resulting system [see Fig. 1(b)] still 
has one more symmetry to be exploited. If one performs 
a symmetric-antisymmetric linear combination of a even and 
peven crea ti n-annihilation operators acting on sites situated 
at the same distance from the QD [see Fig. 1(b)] and a 
bonding-antibonding linear combination of the QD levels \a) 
and |/3), one obtains the MOs |+) and |— ), which couple to 
symmetric and antisymmetric bands, respectively, as shown in 
Fig. 1(c), through hopping matrix elements 



t+ = V2(t' + t") 
t- = V2(t'-t") 



(5) 
(6) 



It is now clear that, at the 2LSU2 point (where v = 1, 
i.e., t' = t", and therefore t_ = 0), a third band decouples 
from the system, leaving MO |— ) decoupled from any elec- 
tron reservoir, although still capacitively coupled to MO |+). 
This immediately answers some of the questions raised above. 
First, in the limit v — > 1, the charge transport changes from 
double- to single-channel. Therefore, the maximum conduc- 
tance will change from 2Gq to Go- Indeed, in the vicinity 
of the 2LSU2 point, for v < 1.0, narrow peaks with height 
2Gq are expected around certain values of gate potential. The 
closer one gets to the 2LSU2 point, the narrower these peaks 
become, their width vanishing at the 2LSU2 point, where the 
maximum of the conductance will be Go pSEl Second, at the 
2LSU2 state, MO |— ), being decoupled from the charge reser- 
voir, will have zero width, i.e., T_ vanishes (for a definition of 
r_, see below and Appendix A). As the gate potential varies, 
and after the level |— ) has crossed the Fermi energy, it will be 
abruptly charged, changing the ground state of the system (re- 
call that MO | — ) has a many-body interaction with the rest of 
the system), and potentially changing its transport properties 
discontinuously. One should expect that the Kondo peak in 
the LDOS, when the system is in the SU(4) state, should dif- 
fer considerably from the Kondo peak occurring in the 2LSU2 
state, given the difference in the number of charge transport 
channels between the two states. The very abrupt change in 
the slave boson's LDOS (seen in Fig. 14 of reference O is 
probably linked to this change from double- to single-channel 
charge transport. If the change in the Kondo resonance is dis- 
continuous or not (at = 7r/4) is still a point to be further 
investigated. 

For completeness, we present here i? m b written in the new 
basis, assuming that U = U' (in Appendix A, a complete 
expression, for the case U ^ V is presented), 



H, 



U + U' 



mb 



^2 n\an\ s + ^2 n Xa nx s 

\,<T A^A',(7 



U' 



n\ a n\>, 



(7) 



where A = ±, a =f or J,, and although our calculations are 
for U — U', both are indicated, with their respective terms. 
It is also useful to define the broadening of the one-body AOs 



| a, f3), as well as for the MOs |±) (coming from their coupling 
to the leads). For the AOs, we have 

r = r a + i>, (8) 

r Q/(9 = 2n(t' 2 + t" 2 )p (E F ) =2Trt 2 oPo (E F ), (9) 

where the factor 2 in Eq. [9] takes in account the left and right 
leads [see Fig.l (a)], and po(Ep) is the LDOS at the first site 
of the leads (here considered the same for the a and j3 chan- 
nels), and T is the total coupling between the interacting re- 
gion and the leads. For the MOs we have 



irtlpoiEp), 
irt 2 _ Po (E F ). 



(10) 
(11) 



Note that, as required, T = T + + T_ (for more details, see 
Appendix A). 

The MO |— ), in the 2LSU2 regime (y = 1), is an elec- 
tronic dark state that cannot interact directly with the con- 
duction electrons. Such state has been extensively studied 
in different systems, namely, as Dicke effecPKH or bound 
states in the continuum (BICs), 25 26 which can produce Fano 
resonances in the conductance. 27 The original Dicke effect in 
quantum optics^ takes place through spontaneous emission, 
by two closely spaced atoms (of the same species), which emit 
a photon into a common environment. 



in. CHARGE OSCILLATION OF MOLECULAR 
ORBITALS: DMRG AND NRG RESULTS 

A. DMRG and NRG (zero and finite-temperature) results 

Figure 2(a) contrasts DMRG results at the 2LSU2 state 
{(f> = 7r/4) obtained for n± vs. V g and n a $ vs. V g , empha- 
sizing the difference between them: the AOs present a smooth 
charging, as a function of V g , and, obviously, n a (V g ) = 
np (V g ), while n± (V g ) show clear discontinuities and very 
diverse behavior. Note that we plot only up to the p-h sym- 
metric point (Vg/U = —1.5), and the data has to be multiplied 
by 2 to reproduce the total charge in the QD, i.e., we are plot- 
ting charge per spin species, as the system is SU(2) invariant. 
Figure 2(b) shows a comparison of n± vs. V g between NRG 
and DMRG, at zero-temperature. The agreement between the 
two techniques is quite good, as can be easily seen. The hop- 
ping and many-body parameters used are indicated in the fig- 
ure and in the caption, and cf> = ir/4 for both panels (2LSU2 
state). As mentioned above, at the 2LSU2 state, level |— ) is 
completely decoupled from the electron reservoir and will ob- 
viously be abruptly charged for some V g value, and, given that 
it is capacitively coupled to level |+), its occupancy will also 
change abruptly. This is clearly reflected in the results for n± 
in Figs. 2(a) and (b). Note that the results in Fig. 2(a) are for 
a smaller value of U/T than those in Fig. 2(b) (the second is 
deeper into the Kondo regime), although this discontinuity in 
the charge occupancy of the MOs still occurs for both cases. 
Note that the log-discretization of the semi-infinite chain used 
by NRG is the one described in Ref.|29l 
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FIG. 3: Left panels: DMRG results similar to the ones in Fig. 2, but 
now for varying to (0.025, 0.05, 0.106, and 0.175, from top to bot- 
tom). The corresponding values for V are 0.0025, 0.01, 0.045, and 
0.1225. U — 1.0 for all panels. For the sake of clarity of presenta- 
tion, just the n+ results are shown. Note that the value of V gc [gate 
potential for charging (discharging) of |— ) (|+))] is strongly depen- 
dent on r, while V g d (gate potential for discharging (recharging) of 
— ) (|+})] is not (its value is practically fixed at V g /U = —1.0). 
Right panels: Solid (black) lines: conductance obtained through the 
Friedel sum rule applied to the corresponding charges in the left pan- 
els. (Red) Dots in panel (c): NRG results for the conductance. Notice 
the good agreement between the two techniques. The narrow peaks 
in the 3 top right panels, located &\V g /U = —1.0, are briefly dis- 
cussed in the text. OBS.: The bottom right panel should be compared 
to Fig. 4 in Ref. [10] The fact that only one NRG result is shown in 
the right-side panels is explained in the text. 



The questions to be answered here are the following: (i) 
what determines the value of V g for which the charging (dis- 
charging) and discharging (recharging) of level |— ) (|+)) oc- 
curs (indicated in Fig. 2(a) with arrows and denoted V gc and 
Vgd)'? (ii) How this value depends on T? (iii) And finally, 
how robust is this phenomenon? I.e., how far from the strict 
2LSU2 regime can we go (i.e., for <f> < tt/4) and still observe 
this gate potential dependent MO charge oscillation effect? 
Why this charge oscillation effect is important will be appar- 
ent soon (see Fig. 6). 

Figure 3, showing DMRG results for different values of 
U /r, answers the second question. I.e., how does the charg- 
ing of MO level |— ) depends on V (for fixed U)l In it, on the 
left side panels, one sees that V gc depends on V much more 
strongly than V g d, as the value of V g d barely changes with T, 
while Vg C decreases as T increases, until V gc = V g d (a full 
explanation of this effect will be given in section V). To un- 
derscore the fact that for <f> = tt/4 charge transport occurs 
only through MO |+), the right panels in Fig. 3 show the con- 
ductance results obtained by using the Friedel sum rule (solid 
curves) G = Go sin 2 (7m + /2), where the results used for 
n + (V g ) are those obtained by DMRG (left panels). The dot- 
ted curve in panel (c) shows the conductance results obtained 




FIG. 4: Temperature (in units of U) dependent NRG results for 
n± vs. V g . Maximum temperature (right pointing triangles) is 
T = 3.118 x 10~ 2 [minimum temperature is zero (circles)]. The ar- 
row in panel (a) indicates direction of temperature increase. Note that 
the MO V g -dependent charge oscillations have been strongly sup- 
pressed at T = 7.888 x 10~ 3 (stars), a value that is of the same 
order of magnitude as the Kondo temperature of a corresponding 
(same parameters) single impurity problem. This suggests that the 
MO charge oscillations arise from a low-energy effect, possibly re- 
lated to the Kondo effect. Parameters are indicated in the bottom 
panel. 



by NRG (through the calculation of the dressed Green's func- 
tion). The good agreement between both techniques indicates 
that, as expected, the conductance is all through the molecular 
orbital |+). The origin of the sharp peak atV g /U = —1.0 
will be explained in section IV, where the LDECA results for 
conductance are presented. Note that the lack of NRG results 
for the other right-side panels is due to convergence issues.®! 

To help answer the other questions, it is important to have 
an idea of the order of magnitude of the energy scale of the 
processes that determine this gate potential dependent charge 
oscillation. The temperature (in units of U) dependent n± 
vs. Vg NRG results in Fig. 4 indicate that these V s -dependent 
charge oscillations have been strongly suppressed at a tem- 
perature which is a very small fraction of U (see temperature 
values in the caption). Indeed, as the charge transport occurs 
through level |+), as shown by the description of the DMRG 
results presented in Fig. 3, the Kondo temperature of the sys- 
tem is in reality that of a single QD with the same T and U 
as that of level |+). This Kondo temperature, calculated us- 
ing NRG, is of the same order of magnitude as the tempera- 
ture for which the V g -dependent charge oscillations have been 
strongly suppressed (T w 7.888 x 1CP 3 , * curve). This in- 
dicates that, as initially thought, these charge oscillations are 
governed by a many-body process related to the Kondo effect, 
or some other many-body effect with a similar characteristic 
energy (as shown in Appendix B, the non-interacting model, 
U' = U = 0, does not show these Vg -dependent charge oscil- 
lations). 
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FIG. 5: (a) n+ vs. V g DMRG results (U = 1.0, V = 0.045, and 
</> = 7r/4) comparing n+ (Vg), solid black line, with n+ (V g + U) 
[(green) + symbols] for V g values around — 0.9[/. The perfect agree- 
ment between the two curves, in the region between discontinuities, 
indicates the simple fact that the charging of the decoupled MO level 
|— ) acts as an extra electrostatic potential, moving up by U the effec- 
tive potential seen by MO level | +) . This is schematically shown in 
(b), where the effective potential is plotted against V g and the (red) 
vertical arrow indicates the charging of MO | — ) and the change by U 
of the effective potential seen by MO | +) (the discharging of MO | — ) 
occurs in the other discontinuity), (c) and (d): DMRG results for n± 
vs. V g for three different values of (0.25, 0.24, and 0.23). This 
shows that the V g -dependent charge oscillations seen at (f>/n = 0.25 
survive away from the 2LSU2 point. The charge oscillations are still 
noticeable down to <f)/n ~ 0.2 (not shown). The parameters are 
t = 0.106, T = 0.045, and U = 1.0. The DMRG calculation 
was done for a cluster with 80 sites and 200 states were kept, and for 
4> < 0.25 only the left-right transformation was used [see Fig. 1(b)] 
The black horizontal arrow in panel (c), when applied to the (green) 
+ symbols curve, schematically shows the displacement presented 
in panel (a). 



B. The effective gate potential model (4> = 7r/4) 

From the point of view of level |+), level |— ), when 
charged, acts as an extra source of electrostatic energy, be- 
sides the one provided by the back-gate (Vg). The particular- 
ities here are that, at the 2LSU2 point (0 = 7r/4), level |— ), 
being decoupled from the electron reservoir, can only accom- 




V \U 

g 

FIG. 6: (a) LDECA results (t = 0.15, T = 0.09, U = 1.0, A = 7, 
and (j> = 0.242) for conductance G [(black) circles], charging of 
the MO levels n± (x (green) and + (blue) symbols for n+ and 
respectively), and charging of the AO levels n a ^ (square (red) sym- 
bols). The mere comparison of G with n± and with n a .fj indicates 
that the conductance variation with V g is intimately correlated with 
that of n±, i.e., understanding the Vg -dependent charge oscillations 
of the MOs j±) is crucial to the understanding of the charge transport 
properties of the system. In addition, by comparison with Figs. 5(c) 
and (d), one notices the qualitative agreement between the DMRG 
and LDECA results for n± as a function of V g . (b) LDECA results 
for spin and charge fluctuations for both MOs. x symbols (green) 
and + symbols (blue) indicate (5S Z ) 2 for MO |+) and |— ), respec- 
tively; diamonds (red) and * (cyan) indicate (Sn) 2 for MO |— ) and 
|+), respectively. A plateau in (5S Z ) 2 , accompanied by simultane- 
ous suppression of charge fluctuations, indicates a level participating 
in a Kondo effect. It is interesting to note that, according to these 
results, MO level |— ), although very weakly coupled to the charge 
reservoir, and therefore still undergoing strong Vg -dependent charge 
oscillations (notice that <j) ~ 7r/4), seems to participate in the Kondo 
effect, mainly close to the p-h symmetric point, where its spin fluctu- 
ations are higher than those for MO |+). The inset shows the results 
for spin-spin correlations between the AO levels \a) and |/3) (verti- 
cal triangles). Note that close to the p-h symmetric point they are 
ferromagnetically correlated. The cartoon at the left side of the inset 
explains why (see text and Reference [34). The conductance struc- 
tures indicated by small arrows in panel (a) are discussed in the text. 
Note that the charge fluctuation for MO |+) [* (cyan)] was divided 
by 2 to fit better in panel (b). 



modate an integer number of electrons (7i_ = 0, 1, or 2). 
Taking this into account, one can define an effective gate po- 
tential as V^* = V g + f/n_ . From this point of view, it is then 
easy to understand the effect of the charging and discharging 
of level |— ): it just modifies the effective gate potential seen 
by level | +) . To emphasize that this is indeed the only effect 
in play for the charging of level |+) at the 2LSU2 state, we 
show in Fig. 5(a) the curves for n + (V g ) and n + (V g + U) for 
Vg/U values around —0.9 [for U/T = 22.0, corresponding 
to Fig. 3(c)]. The horizontal arrow in Fig. 5(c) schematically 
shows the 'superposition' of the two curves shown in Fig. 5(a). 
The matching of the two curves in Fig. 5(a), for n + (V g ) [solid 
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(black)] and n+(V g + U) [(green) + symbols], in the region 
between discontinuities, is perfect, confirming the effective 
potential idea described above: the discontinuous charging of 
the MO |— ) changes by U the gate potential seen by MO |+), 
discharging it by 1 electron (as the total charge in the sys- 
tem is kept constant), and, more importantly, as V g decreases 
further below V gc , the recharging of MO |+) occurs exactly 
the same way [see Fig. 5(a)] as when MO |— ) was empty, at 
V g = V gc + U. As a reference, the effective gate potential V g * , 
as a function of V g , is shown in Fig. 5(b). The question still 
to be answered is what physical process determines V gc (T), 
i.e., the dependence of the value of V gc with the coupling to 
the leads T (as depicted in the left-side panels of Fig. 3). Be- 
fore answering this question, let us concentrate on the other 
two pending issues: how robust is this gate-dependent charge 
oscillation away from the exact 2LSU2 point, and what is its 
effect over the charge transport, i.e., the conductance of the 
system. 



IV. DMRG AND LDECA RESULTS AWAY FROM THE 
2LSU2 POINT 

Figures 5(c) and (d) show DMRG n± vs. V g results for 
three different values of 4>/ir (0.25, 0.24, and 0.23). As ex- 
pected, the discontinuities are smoothened out, but the charge 
oscillations survive. Additional calculations (not shown) in- 
dicate that the oscillations are still noticeable for c/)/tt rj 0.2 
{f It' k, 0.73). Once < 7r/4, the transport problem be- 
comes a double channel one, making its solution through 
NRG much more costly in computational terms. To circum- 
vent this, here we use the LDECA method 16 to present results 
for <fi < 7r/4. Fig. 6(a) shows results for conductance G (* 
(black) symbols), and n± vs. V g for 4>/n — 0.242 (x (green) 
and + (blue) symbols for n + and n._, respectively). The (red) 
square symbols show the charging of the AO levels \a, 0). 
A comparison of the conductance curve with the curves for 
the charging of the MO levels and the AO levels indicates that 
the charging of the AO levels provides very little information 
about the sharp features observed in the conductance. These 
features are obviously more correlated to the charging of the 
MOs. As tj> ps 7r/4, we are still too close to the single-channel 
regime, and the conductance is still too close to those obtained 
in Fig. 3 (right-side panels). The charging results for the MO 
levels are in good qualitative (and even quantitative) agree- 
ment with those for DMRG [as shown in Figs. 5(c) and (d)]. 

It is interesting to contrast the very sharp peak in the con- 
ductance, around the point where level | — ) is discharged, with 
the ones shown in the right-side panels in Fig. 3. In Fig. 5(a), 
the peak reaches one quantum of conductance (Go), while in 
Fig. 3 it reaches just m Go/2). Although the convergence 
of LDECA results as — > ir/4 becomes increasingly dif- 
ficult, the peak at V g = —U, for <j> < 7r/4, has always a 
maximum of Go- We speculate that the difference with the 
(\> = ir/4 results (Fig. 3) stems from the number of chan- 
nels. Although r_ is very small when s» ir/4, there are 
two channels available for charge transport (contrary to the 
situation at <f> = 7r/4). Notice in Fig. 5(a) that the discharg- 



ing of MO |— ) (at Vg/U ~ — 1) occurs over a window of V g 
(with width T_), indicating that level |— ) is in resonance with 
the band, therefore participating in charge transport. LDECA 
results for smaller values of <f> (down to <fr/ir ~ 0.2) indi- 
cate a smooth evolution of the conductance, until values of 
G start to show features above Go, once the double chan- 
nel transport takes hold (and the conductance reaches 2Go for 
some narrow regions of V g ).^M In addition, a comparison of 
the right side tail of the narrow peak, with the corresponding 
right tail of the broad plateau (approximately centered around 
Vg/U = — U /2) shows that the effective potential idea is still 
helpful in interpreting the results (at least for this small de- 
viation from 4> = 7r/4). Note that rt_ (+ (blue) symbols) is 
still very approximately quantized. Finally, it is interesting 
to point out that the small upticks in the conductance, before 
the sharp drop, and after the sharp rise (indicated with arrows) 
have been seen in functional Renormalization Group (fRG) 
calculations, for a similar model, by C. Karrasch et al. [see 
their Fig. 5(a)] P 

In Fig. 6(b), results for spin and charge fluctuations for 
the MOs are shown. As expected, the rightmost conductance 
plateau in panel (a) is associated to a plateau in the spin fluc- 
tuations of MO level |+) (x (green) symbols), and a suppres- 
sion of its charge fluctuations (* (cyan) symbols), which are 
fingerprints of the Kondo effect. Taking in account the com- 
ment above, that the effective gate potential model is still use- 
ful, the narrow peak around V g /U = —1.0 can then be in- 
terpreted as a Kondo state that was 'skipped' due to the sec- 
ond charge oscillation: indeed, when MO |— ) discharges ap- 
proximately 1 electron, the effective potential seen by MO |+) 
abruptly decreases by ~ U, moving the system from the high 
Vg side of the Intermediate Valence regime to its low V g side, 
completely skipping the Kondo regime (note the enhancement 
of charge fluctuations for MO |+) around V g /U = —1.0, and 
a much smaller value for its spin fluctuations, as compared to 
the right-side maximum, indicating that the charge transport 
process here is more akin to sequential tunneling}21 

Around the p-h symmetric point, a Kondo effect with pecu- 
liar characteristics arises, where both MOs participate. This 
can be confirmed by checking the spin fluctuations of MOs 
|+) and |— ). In reality, the spin fluctuations of MO |— ) (+ 
(blue) symbols) are (surprisingly!) more enhanced than those 
of MO |+). Given the capacitive coupling between the two 
MOs, the participation of level |+) in a Kondo effect drives 
level |— ) into also becoming an active degree of freedom in 
the Kondo effect. A somewhat different view of this Kondo 
effect can be taken if one returns to the AO basis and cal- 
culates the spin correlations between AOs \a) and |/3). This 
is shown in the inset in Fig. 6(b): the two AO states, close 
to the p-h symmetric point develop a ferromagnetic correla- 
tion. This can be readily understood by referring to Fig. 1(b): 
instead of performing the symmetric/antisymmetric and the 
bonding/anti-bonding transformations that lead to Fig. 1(c), 
one should skip the second transformation (keeping the AOs). 
This leads to (for v = 1) the schematics shown at the left side 
of the inset in Fig. 6(b), which is well known to lead to ferro- 
magnetic correlations between the AOs, as remarked by Mar- 
tins et a/.^and confirmed in other publications thereafter!^ 
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FIG. 7: (a) Comparison of T (V 9 ) (solid (black) curve) with T (V g + 
U) (dotted (red) curve) for a large window of V g . Notice the cross- 
ing point for the two curves (indicated by an arrow). As can clearly 
be seen, this is the point for which the gain in energy (To) by 'mov- 
ing up' the effective gate potential of the system by U, and going 
back to the Intermediate Valence regime, is larger than the gain in 
energy provided by 'staying' in the Kondo state. It can be clearly 
seen in panels (b) and (c), which show the effective gate potential 
and corresponding results for n±(V g ), respectively, that the crossing 
indicated in panel (a) exactly agrees with V gc , i.e., the point where 
n_ is charged and n+ is discharged. In addition, the second charge 
oscillation (at V g d/U = —1) occurs at the point where To(V g + U) 
crosses To(V g ) for the second time, but now passing under it. In 
this case, it is energetically favorable for the system to 'move down' 
the effective potential, i.e., when the system starts to enter the Kondo 
regime again, it is energetically advantageous for the system to 'skip' 
the whole Kondo regime and go straight to the second Intermediate 
Valence regime. This is accomplished by MO level |— } being dis- 
charged and level |+) being recharged. Parameters used in panels 
(a) to (c) are the same as the ones used in Fig. 3(a). Panels (d) to 
(f) show similar results as in panel (a), but for progressively larger 
values of V (the same parameters as the ones used in panels (b) to 
(d) in Fig. 3). For all the cases analyzed, the crossing of To(V g ) and 
T (V g + U) agrees with V gc . 



V. COMPETITION BETWEEN KONDO AND 
INTERMEDIATE VALENCE EFFECTS (</> = tt/4) 

The origin of the dependence of V gc with the coupling to the 
leads (r), as shown in the left panels of Fig. 3, still has to be 
understood. In Fig. 4 we have shown that the Vg-dependent 
charge oscillations are related to a low energy scale phe- 
nomenon. Another important point that has to be remarked, 
is that the very dependence of V gc with T clearly indicates 
that the MO level |— ) does not get charged the moment it 
crosses the Fermi energy, as V g varies. This conclusion is ob- 
vious from the fact that the position of the Fermi energy is 
the same for all left-side panels in Fig. 3, and still, the level 
|— ) gets charged at different gate potential values. This is a 
consequence of MO level |— ) being strongly correlated to the 
whole system, including the electron reservoir, through MO 
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FIG. 8: Comparison of the DMRG finite size scaling calculation of 
To (see text for its definition), with the Haldane expression for Tk, 
Eq. (12). Curves for TV = 10, 20, and 30, were translated to match 
the TV = 40 curve at V g /U = —0.5. The same translation was ap- 
plied to Haldane's curve. It is clear that, in the region of V g where the 
Haldane expression applies (before it diverges), the DMRG results 
for larger clusters gradually approach Haldane's results. Parameters 
are the same as the ones used in Fig. 3(a). 



|+). Bearing in mind now that the charging of level |— ) at 
a certain V g value results in 'moving back' level |+) from 
a Kondo state to an Intermediate Valence state by increasing 
its effective gate potential by U, one may assume that it is 
energetically more favorable for the system at that V g value 
to have MO level | +) in the Intermediate Valence state than 
in the Kondo state. By using DMRG to calculate the ground 
state energy of a Single Impurity Anderson Model (SIAM) as a 
function of V g , this intuitive idea can be quantitatively probed. 
To achieve that, one just has to remove from the SIAM ground 
state energy, calculated with DMRG, its single-electron con- 
tributions, i.e., -ffband, the gate potential energy, and the local 
correlation energy U (nj, iVf), thus defining the following en- 
ergy, here denoted as T (V g ): 



T (V g ) = E (T = 0, - E Q (T,V g ), 



(12) 



where Eq(T, V g ) indicates the ground state energy (calculated 
with DMRG) for a SIAM with parameters (7 = 1, and hop- 
ping to the band equal to t + = 2<o [for to = 0.025, cor- 
responding to Fig. 3(a)], while Eo(T = 0, V g ) is the corre- 
sponding ground state energy for the same model and same 
parameters, but now with the impurity disconnected from the 
band. Note that, as To (V g ) > 0, it represents the gain in 
energy, by the system, when the many-body correlations are 
extended into the leads (be it by the formation of the Kondo 
state, or any other strongly correlated ground state). 36 We find 
out, as shown in Fig. 7, that this energy To (V^) explains the 
dependence of V gc with T. 

In Fig. 7(a), the energy To 1S plotted as a function of V g , 
and, as done in Fig. 5(a), we compare Tq{V 9 ) (solid (black) 
curve) with T (V g + U) (dot-dashed (red) curve). The 'hor- 
izontal displacement' of Tq by U (dot-dashed (red) curve) 
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mimics the effect over MO |+) (the single impurity in the 
SIAM) of the charging of MO |— ). The arrow in Fig.7(a) 
indicates the crossing point of the two curves, i.e., the point 
for which the gain in many-body energy To, by 'moving up' 
the gate potential of the system by U, and going back to the In- 
termediate Valence regime (dot-dashed (red) curve), is larger 
than the gain in energy provided by 'staying' in the Kondo 
state [(black) solid curve]. What affords the double-MO sys- 
tem (with MO levels |+) and |— )) the ability to put in direct 
competition two states relatively distant in the phase diagram 
(the center of the CB valley and the initial ramp up of the CB 
peak) is the presence of the zero-width strongly correlated MO 
level |— ) that can be charged abruptly and can only contain an 
integer number of electrons (for <j) = ir/4, or approximately 
integer for <f> < 7r/4). Panels (b) and (c) in Fig. 7 show the ef- 
fective gate potential (as a function of V g ) and the correspond- 
ing results for n±(V g ), respectively, for the same parameters 
used to calculate T (U = 1.0, t = 0.025) for the SIAM 
in panel (a), but for the model containing both MO states. 
The crossing indicated in panel (a) (vertical arrow) agrees ex- 
actly with Vg C , i.e., the point where n_ is charged and n + 
is discharged. In addition, the second charge oscillation (at 
Vgd/U = — 1) occurs at the point where Tq(V 9 + U) crosses 
Tq(V 9 ) for the second time, i.e., when MO level |+) starts to 
enter the Kondo effect again (dot-dashed (red) curve in panel 
(a), around v g /U = —1), it is energetically advantageous for 
it to 'move down' (by decreasing the effective potential by 
U, as n_ changes from 1 to 0), thus 'skipping' the whole 
Kondo regime and going straight to the second Intermediate 
Valence regime. This is accomplished by MO level |— ) being 
discharged and level |+) being recharged (and this 'skipping' 
of the Kondo regime results in the narrow peak seen in the 
conductance in Fig. 5(a), centered around V g /U = —1.0). 
Panels (d) to (f) show similar results as in panel (a), but for 
progressively larger values of T (the same parameters as the 
ones used in panels (b) to (d) in Fig. 3). For all the cases an- 
alyzed, the first crossing of To(V g ) and Ta(V g + U) agrees 
with V gc , and, as can be clearly seen, the second crossing oc- 
curs always at the same gate potential V g /U = —1.0 [note 
that in panel (f), the two crossings now coincide, resulting in 
the peculiar charging pattern in Fig. 3(d)]. 

As initially hinted by the finite temperature NRG results in 
Fig. 4, the gate potential dependent charge oscillation effect is 
a low-energy many-body effect; in reality, an interesting com- 
petition between two different many-body states, Kondo and 
Intermediate Valence. The ability of the DMRG calculated 
Tq(V 9 ) to explain this subtle effect motivates us to compare 
To with the well known Haldane expression for Tr-J 2 ° I 37 I 

w ~Kra^) 1/v * |, '•' ll ™■ (13) 

Fig. 8 does just that for a system consisting of a single QD. 
A finite size scaling of To, for DMRG clusters of increasing 
size, N = 10, 20, 30, and 40 (where N is the total number of 
sites in a DMRG one-dimensional cluster containing the im- 
purity plus N — 1 non-interacting sites), shows that T , in the 
V g region around the p-h symmetric point, where the Haldane 



expression applies, is tending to Haldane's Tr-. The param- 
eter values are indicated in the figure. A few points have to 
be stressed here. First, note that the objective of Fig. 8 was 
not to attempt a simulation of the absolute value of Tk (all 
the curves have been translated in energy so that they coin- 
cide with the N = 40 value at the p-h point). The objective at 
this point is to reproduce the dependence of Tr- with gate po- 
tential V g . Second, note that the Haldane expression is by no 
means exact, and is supposed to be applicable only when real 
charge fluctuations of the impurity can be neglected, therefore 
we are not expecting a complete agreement of DMRG results 
with Haldane's Tk- With these two points in mind, and taking 
into account that a larger DMRG computational effort (not at- 
tempted in this work) can greatly improve the accuracy of the 
T results, one can, after a careful analysis of different defi- 
nitions for T), propose a DMRG calculation of To as a useful 
guide to experimentalists for all the relevant regimes of the 
single impurity Anderson modelPS 



VI. CONCLUSIONS 

In this work, motivated by recent interest in CNT QD 
charge transport measurements, the authors use three distinct 
numerical techniques to study the 2LSU2 regime and present 
details of its properties, which hinge in two important aspects: 
first, a transformation from atomic to molecular orbitals is 
needed to understand the numerical results for charge trans- 
port. Second, this transformation unveils an interesting gate- 
potential-dependent charge oscillation involving two strongly 
interacting molecular orbitals, one of which, MO |+), is cou- 
pled to the band, while the other, MO |— ), is totally decou- 
pled from it (at the 2LSU2 state), while capacitively coupled 
to | + ) . In doing that, a series of unanswered questions asso- 
ciated to the 2LSU2 regime are clarified. First among them 
is the fact that at the 2LSU2 regime the maximum conduc- 
tance passes from 2Go to Go, as the charge transport changes 
from double- to single-channel. In addition, all the charge 
transport properties specifically associated to the 2LSU2 state 
are shown to be robust, in the sense that they can in principle 
be observed even when there is no total mixing of the orbital 
quantum numbers (<p < 7r/4). Finally, by numerically defin- 
ing the quantity T , it is shown that this gate dependent charge 
oscillation results from the competition between the Kondo 
and Intermediate Valence states. These two states, separated 
by a wide crossover region, can be effectively put in contact 
by the presence of a strongly correlated dark state (i.e., a zero 
width state), the MO |— ). By comparing our numerical re- 
sults for T with Haldane's T K , we speculate that the DMRG 
results obtained by this relatively easy computational calcula- 
tion may be of interest to experimentalists in the analysis of 
other systems. A more refined analysis of To is now being 
pursued and will be presented elsewhere. In addition, it will 
be the subject of a forthcoming work to understand how the 
phenomena observed in the 2LSU2 regime change for U' ^ U 
and when the degeneracy of the AOs \a) and |/3) is lifted. 
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decouples from the system. The Hamiltonian i?hyb in this 
new basis takes the simple form, 



hyb 



(A6) 



A;A';ct 



Note that the couplings are now renormalized by a factor \[2 
(the factor 2 accounts for double-counting). It is easy to show 
that the band Hamiltonian -ffband nas i ts li near chain struc- 
ture preserved under this transformation, as crossing terms 
between different chains cancel out. The system in this new 
basis is depicted in Fig.[TJb). 



Appendix A: Canonical transformations 

The system studied in this work consists of four different 
Fermi seas connected to a two-level CNT QD, schematically 
shown in Fig. [TJa). Canonical transformations, which elim- 
inate part of the complexity of the problem, 39 40 can be per- 
formed over the system's Hamiltonian, given by Eqs.|2]to|4] 
In this Appendix, we briefly discuss these canonical transfor- 
mations and how to calculate the Green's functions and the 
conductance. 



1. Left-Right symmetry 

The first transformation exploits the left-right symmetry. 
As a consequence of this symmetry (see Eq. [4j, just the sym- 
metric combination of left-right creation and annihilation op- 
erators in the a (/?) band are coupled to the a ((3) orbital of 
the CNT QD. This can be easily verified by simply rewriting 
Eq(4]as 



hyb 



A:A' 



CX'la + CA'-lcr) + H.C. 



(Al) 



where A and A' can take values a or /3. 

Now, constructing the symmetric and anti-symmetric com- 
binations of the c\±i a operators (where A = a or /3) for all 
i < 1, we can define a new basis composed of the following 
orthonormal states, 



c j3ia 
C /3icr 



1 

71 

1 

71 
i 

71 
i 

71 



[C/3i. 



C/3- 



(A2) 
(A3) 
(A4) 
(A5) 



Note that this change affects only the leads' fermion opera- 
tors, keeping the CNT QD operators unchanged. In this new 
basis, the symmetric (even) combination will remain coupled 
to the CNT QD, while the anti-symmetric (odd) combination 



2. Using a- ft Symmetry 

As the CNT QD sits at the center of the system (equiva- 
lent to site 0) it is not affected by the previous transformation. 
Now we can take advantage of the a-(3 symmetry to simplify 
even more the problem. Note that under this new a-[3 trans- 
formation, the many-body terms will be also affected. The 
proposed transformation consists of using the symmetric and 
anti-symmetric combinations between states c e aia and Cp ia 
[see Fig. |TJb)], and keeping the decoupled (odd) states un- 
touched. The new basis for the even subspace will be (where 
now i > 1, i.e., i takes only positive values) 



\_ C aia C l3ia\ 



he 



1 

71 

— [data + dj3cr] , 
1 

71 



[dan — dp a ] 



(A7) 
(A8) 
(A9) 
(A10) 



where the new operators d +a and act on the a, f3 levels of 
the CNT QD. States associated to operators d +a and d^ a are 
referred to as molecular orbitals (MO) in the main text. As in 
the previous transformation, the band Hamiltonian iTband re- 
tains its independent double-linear-chain structure in this case 
too. It is important to note that the total number of particles in 
the CNT QD is conserved under this transformation. For the 
special case U 1 = U, the many-body terms in the Hamiltonian 
-ffmb just depend on the total number of particles. Therefore, 
its form will not change under this transformation: 



U 



(All) 



cr:A=± 



where N T = J2a( n +° + n -<y)- 

The Hamiltonian describing the connection between the 
CNT QD and the leads takes the form, 



hyb = E [ f + d +<r a i° + l - d- a bl<j + H.c. 



, (A12) 



where t± is given by Eqs.|5]and|6] An schematic of the effect 
of this transformation is shown in Fig.[TJc). 
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It is therefore easy to see that in the fully symmetric case 
(*" = t'), and using that t_ = V2(t' - t"), the MO |-) 
will be only capacitively coupled to the system, as i_ = 0. 
As mentioned in the main text, MO |— ) is an electronic dark 
state. This effect is also known in the literature as Dicke effect 
or bound state in the continuous (BIC). In the absence of in- 
teractions (U' = U = 0) the MO |— ) is completely decoupled 
from the charge reservoirs and appears as a delta resonance in 
the LDOS of the CNT QD, moving with V g . For U' = U ^ 0, 
this state couples capacitively to MO |+) and contributes to an 
effective gate potential (acting on level |+)) given by 



v; = v g 



Un. 



(A13) 



as explained in the main section of the manuscript. 

For the sake of completeness, the many-body Hamiltoninan 
for the case U' ^ U is provided in the |±) basis 



Hi = 
H 2 = 



Hi + H2 + H3 + H4 



H 3 
Hi 



U' 


+ 


U 




2 




V 




u 




2 




U' 




u 


4 



(A14) 

{nxt n H + n ^n- X i) , (A15) 



A=+, 



E (4 t 4 



A=+, 



E 

A=+- 

U' ( n+f n_t + n+i n_j.) 



(A16) 

(A17) 
(A18) 



Note that, in the case U' — U, and H3 vanish, and we re- 
cover Eq. All In addition, in the general case ([/' 7^ U), 
MOs |±) will not only be capacitively connected but also 
spin-spin correlated, with the character of the correlations 
(ferro or antiferro) being dependent on the relative values of 
U and U'. This points to a possibly very rich phase diagram, 
which will be fully studied in a forthcoming publicationPS 



3. Green's functions and conductance 

In the main text, we studied this system using three differ- 
ent numerical techniques that require different formalisms to 
calculate the transport properties. We start by explaining how 
to calculate the conductance using LDECA. 

In general, transport properties are related to Green's func- 
tions that propagate electrons from one lead to another (left 
to right, for example). For the system described in Fig. 1, the 
conductance can be written as, 



G = 4vr 2 |. 9o | 2 t 4 [ Po (E F )f 



(A19) 



where 



|ffo| 2 = \G a -i, a i\ 2 + \Gi3-iJi\ 2 

+ |G (3 _ 1 , Ql | 2 + |G Q _ 1 ^ 1 | 2 , (A20) 

and Ga-i.a'i stands for the Green function evaluated at the 
Fermi level Ep, 



where we keep the notation used in Eq. A2 for the cre- 
ation/annihilation operators and the double bracket ((; )) is the 
standard notation for Green's functions. Note that for simplic- 
ity we removed the spin subindex from the later equations. As 
our Hamiltonian does not have spin-flip terms, all terms in the 
propagators in this section have the same spin subindex and 
summation over repeated a indexes is assumed. 

Now we want to rewrite Eq. |A19|using the canonical trans- 



formation given by Eqs. A2 to |A5| By replacing in Eq. A21 
the anti-transformation, we obtain two types of propagators, 



A-1.A1 = ^{(( CA1 ; CA1 t »-((^ 1 ;c^ t ))},(A22) 

(A23) 



Ga-i,a<i = ~m,v&)) 



LDECA is a numerical technique where Green's functions 
can be easily evaluated. We solve, using LDECA, the system 
represented in Fig. |TJb). Then, using Eqs. A19 A20 A22 and 
|A23| we can calculate the conductance as shown in Fig. [6J 

We use two numerical techniques more, NRG and DMRG, 
but for the transport properties we limit these techniques to 
analyzing the case t" — t' (cf> — 7r/4). For this special 
case, we use all the symmetries of the system [see Fig.[TJc)]. 
Within NRG, we calculate the Green's function propagator at 
the MO |+). Then, using the equations of motion for Green's 
functions,^ we left-right anti-transformation, obtaining the 
necessary propagators to evaluate Eqs. A19 and A20 For 



the DMRG technique we use a different approach. After the 
canonical transformations, using the left-right and a- (5 sym- 
metries, we have that just MO |+) is connected to a reservoir. 
Using the fact that the system is a Fermi liquid, we can apply 
the Friedel sum rule, i.e., using the charge at the MO level 
|+), the conductance can be calculated by the equation 41 



G = 2G 



(A24) 



This procedure cannot be used if <fi ^ ir/4, the conduc- 
tance cannot be calculated as just described. However, the 
occupancy of the CNT qD levels a and /3, as well as MO |+) 
and |— ), can still be evaluated using DMRG. 



Appendix B: Non interacting case 

In the absence of interactions (U = 0), the physical prop- 
erties of the system described in Fig. 1, as for instance, the 
local density of states (LDOS) and the conductance, can be 
calculated exactly. 

We can start by analyzing the LDOS at each level of the 
QD. As we mentioned in Appendix [A] a canonical transforma- 
tion can be performed mixing levels a and /? of the QD into 
two new levels, a symmetric MO |+) and an anti-symmetric 
MO I -) . The LDOS at any of the AO levels of the QD can be 
calculated as the superposition of the LDOS at MO |+) and 
MO |-) (since p a = pp), 



Ga-i.a'i = «CA'i;c] k _ 1 ))(a; = E F ), 



(A21) 



PaH 



(Bl) 
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Vg/U 



0.2 -0.1 0.1 0.2 0.3 
Vg/U 



FIG. 9: (a) LDOS of levels a and /3 for V s = 0, to = 0.15, and 
different values of <f>, When <j> — > n/4, MO |— ) gradually decou- 
ples from the band and a sharp resonant state develops at V g in the 
LDOS. (b) LDOS for MOs |±) for = 0.24tt. (c) Conductance as 
a function of V g calculated for to = 0.15 and different values of <f>. 
When cj> = 7r/4 the MO |— ) is totally decoupled and the conduc- 
tance reduces by half, indicating that just one channel is active, (d) 
Charging of levels a, /3, and MOs | ±) as a function of V g for the case 
4> — 0.2tt. As MO |— ) is weakly coupled, its charging resembles a 
step function. 



with A = a or j3. 

For U = V = 0, MOs |+) and |— ) are totally independent. 
Their LDOS is associated to the corresponding coupling with 
the contacts. Each of these levels will have a width given by 



T s (uj) = 7T^ 2 p (w), 



(B2) 



where 6 stands for |+) or | — ) and po is the LDOS at the leads. 
As shown above, in the main text and in Appendix A, when 
— >• 7r/4, the anti-symmetric state gradually decouples from 
the charge reservoir (as t- — >• 0) and its LDOS becomes a 
delta function. When = 7r/4 the MO |— ) is completely 
decoupled and p a = pp is a superposition of a delta function 
and a resonant peak with width IV. A similar effect can be 
observed in other DQD systems^!!] or j n a double Rashba 



nn: 



,24 



Figure [9| a) shows p a and pp for different values of 0. The 
LDOS was calculated using the Green's function equations 



of motion formalism, which is exact in the limit of U = 0. 
Fig.gb) shows the LDOS for the MOs |+) and |-) for = 
0.247T and V g = 0.0. For this value of 0, close to it/ 4, the MO 
|— ) is almost decoupled and therefore appears as a narrow 
resonance located at V g . 

The formula for the conductance is given by^ 



G(V g 



2G 



io[t 2 V*(l + sin 2 20) + t cos 4 20] 
T+(V g ) T_(V g ) 



(B3) 



T±(V g ) = f (l ± sm2cf>) 2 +t 2 V g 2 
where Go = 2e 2 /h and to — 4tg. 

When levels a and f3 are not mixed (t" 



(B4) 

0, i.e., cj) = 0), we 
can see that the maximum of the conductance occurs for V g — 
0, taking the value G = 2Go, indicating that both channels 
are active. 

As seen above, for = 7r/4, one of the channels is decou- 
pled and the maximum of the conductance is Go- For values 
of <j> close to 7r/4 Eq B4 can be written as a superposition of 



two Breit-Wigner line shapes 



G(0 ~ vr/4) ~ G 



26 27 



where 



A+ 
A_ 



At 2 /t, 

2tl/t [1 - sin 2</>] 



(B5) 



(B6) 
(B7) 



Using Dicke effect terminology, the first term in Eq. |B5] cor- 
responds to the contribution of the super-radiant state (sym- 
metric channel |+)) with a width A+, and the second term 
corresponds to the sub-radiant state (anti-symmetric channel 
|— )) with a width A_. 

In Fig. |9]c) we show the conductance for different values 
of 0. For = both channels are active and have the same 
width. By increasing 0, the MO |— ) starts to decouple and 
this is reflected in a narrower line in the conductance super- 
posed on the wider line given by MO |+). When = 7r/4 
the MO |— ) is completely decoupled and does not give any 
contribution to the conductance, whose maximum is reduced 
from 2 Go to Go. 

Panel (d) of Fig. [9] shows the charging for the case = 
0.27T. Note that the charging of MO |— ) is almost a step func- 
tion. When = tt/4, the charging of levels a and (3 will occur 
discontinuously. 
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